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We report the first detailed comparisons of tlie rates and spectra of neutral-current neutrino 
interactions at two widely separated locations. A depletion in the rate at the far site would indicate 
mixing between v^ and a sterile particle. No anomalous depletion in the reconstructed energy 
spectrum is observed. Assuming oscillations occur at a single mass-squared splitting, a fit to the 
neutral- and charged-current energy spectra limits the fraction of f^ oscillating to a sterile neutrino 
to be below 0.68 at 90% confidence level. A less stringent limit due to a possible contribution to the 
measured neutral-current event rate at the far site from v^ appearance at the current experimental 
limit is also presented. 

PACS numbers: 14.60.St, 14.60.Pq 



Several experiments observing charged-current inter- 
actions of neutrinos have provided compelling evidence 
for i^^ and i/g disappearance as the neutrinos propagate 
from the point of production [l|, 0, 0, H, Q ■ The Super- 
Kamiokande experiment has reported extensively on the 
disappearance of i'^ produced in the atmosphere [2]. 
Measurements of solar i^e showed that the disappear- 
ance of those neutrinos is due to matter enhanced conver- 
sions [3|. The KamLAND reactor experiment provided 
clear evidence for i7g mixing [J| . 

These results are conventionally interpreted as mix- 
ing among the active neutrino flavors that couple to the 
electroweak current. Precise measurements of the Z bo- 
son decay width indicate there are only three light active 
neutrinos [6|], but they do not exclude the existence of 
"sterile" neutrinos, z^s, that do not couple to the elec- 
troweak current. Sterile neutrinos could help resolve sev- 
eral outstanding problems in particle physics and astro- 
physics. For example, sterile neutrinos with masses at 
the eV energy scale can participate in the seesaw mech- 
anism to introduce neutrino masses f7| and can also aid 
in heavy element nucleosynthesis in supernovac 8|. The 
SNO experiment has shown that the total flux of active 
neutrinos from the Sun agrees with the expectation from 
solar models [9] , thereby limiting the extent to which the 
first or second neutrino mass eigenstates could couple to 
a sterile neutrino. While the Super-Kamiokande exper- 



iment excludes pure f^ 



and favors pure Vf^ 



oscillations in its analysis of atmospheric neutrinos, an 
admixture of the two possibilities is allowed ^lOJ and has 
attracted considerable attention in the literature [11| . 

The MINOS experiment has reported a significant 
deficit of i^^ at its far detector relative to the near de- 
tector through measurement of the rate of v^ charged- 
current (CC) interactions [^, ll4|. If this deficit is due 
solely to conversions of i^^ to fr + i^e, then the rate of 
neutral-current (NC) interactions at the far detector re- 
mains unchanged from the non-oscillation prediction. Al- 
ternatively, if any ly^ convert to a sterile state, then the 
NC rate would be suppressed and the reconstructed en- 
ergy spectrum would be distorted. In this Letter we re- 
port the first measurement of the total active neutrino 
rate using a precisely known long baseline and neutrinos 



'Deceased. 



produced with an accelerator. The reconstructed energy 
spectra for NC and CC interactions are used to limit the 
fraction of v^ converting to Vs by fitting them to a model 
of oscillations between j/^ , i/,- , i/g , and Vg dominated by 
the atmospheric mass-squared splitting. 

The neutrino beam is produced using 120 GeV/c pro- 
tons from the Fermilab Main Injector incident on a 
graphite target, which is followed by two magnetic fo- 
cussing horns. The neutrino energy spectrum can be 
changed by adjusting the horn current or the position of 
the target relative to the horns. The flavor composition 
of the beam is 92.9% v^, 5.8% 17^, and 1.3% Vg + Ve- In 
this analysis the v and V are assumed to oscillate with the 
same parameters. The data used in this analysis come 
from the low energy beam configuration whose peak neu- 
trino energy is 3.3 GeV 0, [iji with an exposure of the 
far detector to 2.46 x 10^° protons on target. 

The MINOS near detector is located 1.04 km down- 
stream of the target, has a mass of 0.98 kt, and lies 103 m 
underground at Fermilab. The far detector is 734 km 
downstream of the near detector, has a mass of 5.4 kt, 
and is located in the Soudan Underground Laboratory in 
Minnesota, 705 m below the surface. The fiducial masses 
used for the near and far detectors are 27 t and 3.8 kt 
respectively. 

The MINOS detectors are steel scintillator track- 
ing calorimeters [Ij]. The vertically oriented detector 
planes are composed of 2.54 cm thick steel and 1 cm 
thick plastic scintillator. The scintillator layer is com- 
prised of 4.1 cm wide strips with each strip coupled via 
wavelength-shifting fiber to one pixel of a multi-anode 
photo- multiplier tube [ij,[i5'|- The ncar(far) detector is 
magnetized to an average toroidal field of 1.3(1.4) T. 

Hadronic showers resulting from NC interactions gen- 
erate scintillation light in an average of 12 strips for 1 
GeV of deposited energy. Events must have at least 4 
strips with signal in order to be considered in the analy- 
sis. Individual scintillator strips are grouped into either 
reconstructed tracks or showers, which are combined into 
events. The vertex for each event is required to be suffi- 
ciently far from any edge of the detector to ensure that 
the final-state hadronic showers are well contained within 
the fully sampled portion of the detectors. 

The near detector data are used to predict the number 
of expected events in the far detector, but the ability to 
make this prediction is complicated by the high rate envi- 
ronment at the near detector. At an intensity of 2.2 x 10^^ 



protons on target, an average of 16 neutrino interactions 
are produced in the near detector for each spiU ^. The 
reconstruction program separates individual neutrino in- 
teractions that occur within the same spill. This initial 
pass overestimates the number of NC interactions having 
reconstructed energy, E'roco, < 1 GeV by 36%. Additional 
selections making use of event topology and timing are 
then used to decrease this background. Events must be 
separated by at least 40 ns, and events that occur within 
120 ns of each other must have verticies separated by at 
least 1 m in the longitudinal direction [l6|. After apply- 
ing these criteria, the remaining background from poorly 
reconstructed events with Sroco < 1 GeV is 7%. 

The rate of neutrino interactions from the neutrino 
beam in the far detector is much lower than in the near 
detector, with approximately 1 interaction for every 10 
spills. Interactions from the beam neutrinos are identi- 
fied using a window around the GPS time stamp of the 
spills of — 2 /IS < i < 12 /xs where i = is the expected 
start time at the far detector of the 10 /is spill. Given 
the low rate of neutrino interactions in the far detector, 
spurious events that are coincident with the beam spills 
from noise, cosmic-ray muons, or poor event reconstruc- 
tion can introduce backgrounds to the analysis. Addi- 
tional criteria are used to remove such events, leaving a 
residual background of < 1% of the signal \ll\ . 

Charged-current interactions are identified by the pres- 
ence of a track that may or may not be associated with 
a shower. Neutral-current interactions typically have a 
single hadronic shower, although the reconstruction may 
identify a track in the event; such tracks could come from 
pions, but are mostly reconstruction artifacts. An event 
is classified as NC-like if it has a reconstructed shower, 
is shorter than 60 planes, and has no track extending 
more than 5 consecutive planes beyond the shower [18l |. 
Distributions of these event-topology parameters for near 
detector events are shown in Fig. [TJ The principal back- 
ground in the spectrum of NC-like events comes from 
highly inelastic z^^-CC interactions. The -Ereco spectrum 
of NC-like events in the near detector is shown in Fig. [51 
The distributions in Figs. [1] and [2] show good agreement 
between the data and Monte Carlo simulation. 

The Monte Carlo simulation is used to make an initial 
estimate of the ratio of event yields in the far and near 
detectors as a function of -Eroco- This ratio is multiplied 
by the observed energy spectrum in the near detector to 
produce a far detector prediction of the NC-like event 
spectrum. The true energy of the simulated neutrinos in 
each reconstructed energy bin of the prediction is used 
to determine the effect of oscillations for that range of 
reconstructed energy. To avoid biases, the methods for 
identifying NC-like events and predicting the far detector 
spectrum were developed and tested using only the near 
detector data and Monte Carlo simulation, and the anal- 
ysis procedures were finalized prior to examining data in 
the far detector. 

Figure [3] shows the measured and predicted E^^co spec- 
tra at the far detector. The spectra are compared using 
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FIG. 1: Distributions of event-topology parameters used to 
separate NC-like from CC-like events. Data from the near de- 
tector (solid points) are shown superposed on the total Monte 
Carlo expectation. The hatched distribution shows the //^-CC 
background as determined by the Monte Carlo simulation. 
The systematic uncertainty for the Monte Carlo expectation 
is shown by the shaded band. 



a statistic, R, which expresses the agreement between 
the predicted and observed number of events in the far 
detector: 



R 



Nb. 



Bcc 



•S'nc 



(1) 



where, within a given energy range, A^Data is the mea- 
sured event count, Bcc is the extrapolated CC back- 
ground from all flavors, and S'nc is the extrapolated num- 
ber of NC interactions. The values of S'nc and contribu- 
tions to Bcc are calculated in the framework of three 
neutrino oscillations and are shown in Table IH Because 
the disappearance of Vf^ occurs mainly for true neutrino 
energies < 6 GeV [l2|, the data are separated into two 
samples. Events with -Ercco < 3 GeV are grouped into 
a low-energy sample while events with 3 GeV < Sroco < 
120 GeV are grouped into a high-energy sample. The 
median true neutrino energies of the low and high energy 
samples are 3.1 GeV and 7.9 GeV respectively. The val- 
ues of R calculated for these ranges in i^reco are shown 
in Table HI In the region with i^rcco < 3 GeV, R differs 
from 1 by 1.3cr. Over the full energy range, — 120 GeV, 
the depletion of the total NC event rate is limited to be 
below 17% at 90% confidence level. 

The principal sources of systematic uncertainty in R 
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FIG. 2: The reconstructed energy spectrum for NC-like events 
in the near detector. The data (sohd points) and the Monte 
Carlo expectation including systematic uncertainties (solid 
histogram with shaded band) are shown. The hatched dis- 
tribution shows the expected z/p-CC background. 



TABLE I: Values of A'^oata, •S'nc, and the contributions to Bcc 
for various reconstructed energy ranges. Also shown are the 
values of R. The numbers in parentheses are calculated in- 
cluding Vc appearance at the upper limit discussed in the text. 



Er, 



(GeV) Anata Snc 



Br 



S? 



B"r 



3 
120 



100 
191 



101.1 
98.0 



11.2 
64.2 



1.0 
3.5 



1.8 (9.3) 
11.8 (24.6) 



0-3 i? = 0.85 ±0.10 ±0.07 (0.78 ± 0.10 ± 0.07) 

3-120 7?= 1.14 ±0.14 ±0.10 (1.02 ± 0.14 ± 0.10) 

0-120 71 = 0.99 ±0.09 ±0.07 (0.90 ± 0.09 ± 0.08) 



are listed in Table HTl The absolute scale of the hadronic 
energy is known to within 12%, of which 10% reflects 
uncertainties in the final-state interactions in the nu- 
cleus and 6% results from uncertainty in the detector 
response to single hadrons. The relative calibration of 
the hadronic energy between the two detectors has an 
uncertainty of 3% [5| , and the relative normalization be- 
tween the detectors has an uncertainty of 4%. The un- 
certainty in the near detector event count due to the se- 
lection criteria is 15% for i^reco < 0.5 GeV; 3% for events 
with 0.5 GeV < i?rcco < 1 GeV; and is negligible for 
-Eicco > 1 GeV. The effect of these uncertainties on R is 
shown in Table [TTl 

The uncertainty on the size of the j/^-CC background 
was determined by comparing the near detector NC-like 
reconstructed energy spectrum from the low energy beam 
configuration used in this analysis with the spectra from 
three other beam configurations with higher average neu- 
trino energy. In each reconstructed energy bin, i, of the 
low energy beam the total number of events is the sum 
of the NC and CC interactions, N, = NC, + CC,. The 
quantity rf'^{rf'^) is defined as the ratio of the number 
of NC(CC) interactions in each energy bin in an alter- 
native beam configuration to the corresponding number 
in the low energy beam configuration. The value of CC; 



TABLE II: Sources of systematic uncertainties considered in 
this analysis and their effect on R. 





0-3 GeV 


3 - 120 GeV 


Absolute -Bhad 


± < 0.01 


±0.05 


Relative Shad 


±0.03 


±0.04 


Normalization 


±0.04 


±0.08 


Near detector selection 


±0.02 


- 


Vfi-CC background 


±0.03 


±0.01 


Total: 


±0.07 


±0.10 



can be calculated from the spectrum in another beam, 



,NC 



CC,; = 



Nr - Nf 



,NC 



^cc 



(2) 



where N/^ is the total number of events observed in the 
alternate beam configuration. The values of rf^ and r^*-' 
are taken from the Monte Carlo simulation. The uncer- 
tainty in the i^n-CC background is taken as the difference 
between the uncertainty-weighted average value of CCi 
measured using the different beam configurations and the 
value predicted by the Monte Carlo simulation. That 
difference is consistent within 15% for all reconstructed 
energies. The size of the i^^-CC background at the far 



detector depends on the parameters for Vu 



oscil- 



lations used in the prediction. The MINOS measured 
values of Am§2 = 2.43 x lO^^ eV^/c^ and 6*23 = 7r/4 ^ 
were used for the prediction, and variations within the la 
range of these parameters change the f^-CC background 
in the far detector by less than 10%. 

Because the selection criteria identify i/g-CC interac- 
tions as NC-like with nearly 100% efficiency, the back- 
ground from Pf, inherent in the beam and Vf^ — > Ve os- 
cillations is also considered. An upper limit for the i^e- 
CC rate in the far detector was estimated using the nor- 
mal mass hierarchy with 6112 = 0.61 rad, 6113 ~ 0.21 rad, 
S = Sn/2 rad, Am^, = 7.59 x IQ-^eY^/c^, and Ato§2 = 
2.43 X 10^3 eVVc* EUl- The choice of 6*13 corresponds 
to the 90% confidence level upper limit for the chosen 
Am^2 value [l^. The contribution to Bqc from t/g and 
the values of R in the different energy ranges under these 
assumptions are shown in Table ID 

The data shown in Fig. [3] can be combined with the 
data from CC interactions to determine whether the pre- 
viously observed Vf^ disappearance is due solely to os- 
cillations between the active neutrinos, or if oscillations 
between active and sterile neutrinos also occur. To de- 
termine the fraction of f^ that have converted to a sterile 
state, the data are fit to a model that assumes oscillations 
between v^^Vt^ and Vs occur at a single mass-squared 
splitting. The probabilities for v^ to remain j/^ or con- 
vert to v^ are 






= I- a^sm^ {1. 27 l^m^L/E), and 

= asS\T?{l.2mm'^L/E), (3) 



where Am,^ is the atmospheric mass-squared splitting in 
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FIG. 3: Spectrum of observed NC-like events in the far de- 
tector with predictions for the two oscillation hypotheses de- 
scribed in the text. The filled regions in each bin indicates 
the systematic uncertainty in the predicted rates. 



with x^ = 46.5 for 43 degrees of freedom and fs < 0.68 at 
90% confidence level. The fit includes the systematic un- 
certainties in Table ITTl as nuisance parameters. Including 
electron neutrino appearance at the previously discussed 



upper limit results in fs — 0.43]]Q27(stat.-fsyst.) with 



X 



= 46.6 and /, < 0.80 at 90% confidence level. 



In summary, we have reported the first measurements 
of neutrino neutral-current rates and spectra in an ac- 
celerator long baseline neutrino experiment. The rates 
at the near and far detectors are consistent with expec- 
tations from decay kinematics and geometry, providing 
new support for the interpretation of muon neutrino dis- 
appearance as oscillations among the three active neutri- 
nos. This result provides the best limits to date on the 
fraction of muon neutrinos which may convert to sterile 
neutrinos in oscillations associated with the atmospheric 
mass-squared splitting. 



eV'^/c^, L — 735 km, E is the energy of the neutrino in 
GeV, and a^ and a^ are phenomenological parameters 
related to the mixing angles. A simultaneous fit to the 
Ercco spectrum in Fig.[3]and the i^^— CC energy spectrum 
yields the energy independent fraction of i/^ that oscillate 



to i/. 



P. 



l-Pu 



0.28: 



-0.25 
-0.28 



(stat.-l-syst.), (4) 
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